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Motivated by a recent experiment on nonlinear tunneling in a suspended Carbon nanotube con- 
nected to two normal electrodes [S. Sapmaz, et ai, Phys. Rev. Lett. 96, 26801 (2006)], we 
investigate nonequilibrium vibration-mediated sequential tunneling through a molecular quantum 
dot with two electronic orbitals asymmetrically coupled to two electrodes and strongly interacting 
with an internal vibrational mode, which is itself weakly coupled to a dissipative phonon bath. For 
this purpose, we establish rate equations using a generic quantum Langevin equation approach. 
Based on these equations, we study in detail the current- volt age characteristics and zero-frequency 
shot noise, paying special attention to the advanced or postponed of the appearance of negative dif- 
ferential conductance and super-Poissonian current noise resulting from electron-phonon-coupling 
induced selective unidirectional cascades of single-electron transitions. 

PACS numbers: 85.65.-|-h, 71.38.-k, 73.23.Hk, 73.63.Kv, 03.65.Yz 



I. INTRODUCTION 

Recent progress in nanotechnology has facilitated 
the fabrication of single-electron tunneling devices 
using organic molecules. In turn, this has given 
rise to a large body of experimentalii^i^i^i^i^i and 
theoretical^ii^iiiii^iHiiia^iiaiilii^ work concerning 
phonon-mediated resonant tunneling through a quantum 
dot (QD) with strong coupling to an internal vibra- 
tional (phonon) mode (IVM). In particular. Carbon 
nanotubes (CNT) have recently become the focus of 
much research interest because electronic transport 
measurements show that they demonstrate perfect 
signatures of phonon-mediated tunneling, e.g. stepwise 
structures in the current-voltage characteristics having 
equal widths in voltage and gradual height reduction by 
the Franck-Condon (FC) factor,2iI 

More interestingly, the experimental measurement of 
S. Sapmaz, et al7 has revealed some more subtle trans- 
port features in a suspended CNT connected to two 
electrodes, particularly the appearance of a weak nega- 
tive differential conductance (NDC) at the onset of each 
phonon step followed by a sudden suppression of cur- 
rent (a strong NDC) at a finite bias voltage after several 
steps for a long CNT sample. Theoretically, the reason 
for the weak NDC is quite clear in that it is ascribed 
to the combined effect of strong electron-phonon cou- 
pling (EPC) and low relaxation, i.e. an unequilibrated 
phonon (hot phonon), and strongly asymmetric tunnel- 
couplings to the left and right electrodesi^ii^ii^ Further- 
more, McCarthy, et al. have theoretically predicted the 
catastrophic current decrease, but their calculations are 
based on the presumption that the EPC is dependent 
on the applied bias voltage. It is appropriate to explore 
yet another explanation of the origin of the strong NDC 



without such an assumption. 

Up to now, most theoretical works have focused on the 
case of a single level coupling to the phonon mode. How- 
ever, it is believed that this catastrophic current decrease 
is intimately related to the tunneling processes involving 
two electronic energy levels. Nowack and Wegewijs have 
considered vibration-mediated tunneling through a two- 
level QD with asymmetric couplings to an IVM.^^ Al- 
beit their results display strong NDC through a compe- 
tition between different FC tailored tunneling processes 
of the two levels, their calculations do not fully resolve 
the overall features of the experimental data in Ref. 0- 
Following a suggestion by Hettler, et al.^ that the inter- 
play of strong Coulomb blockade and asymmetric tunnel- 
couplings of two levels to electrodes can lead to a strong 
NDC under certain conditions without EPC, we employ 
a model with two electronic orbitals having asymmetric 
couplings to the electrodes and with both strongly in- 
teracting with an IVM. We establish approximate rate 
equations at high temperature to describe resonant tun- 
neling incorporating the unequilibrated phonon effect. 
Our results are in good qualitative agreement with the 
experiment data of Ref. 0: in the weak bias voltage re- 
gion where the first molecular orbital (MO) is dominant 
in electronic tunneling, the current shows the FC-type 
steplike structure and weak NDC at the onset of each 
step; while with increasing bias voltage the second MO 
becomes dominant in tunneling, inducing a sudden rapid 
reduction of current whether there is EPC or not. Fur- 
thermore, our studies predict that the onset of strong 
NDC as a function of voltage may differ from the value 
at which the second MO becomes active in resonant tun- 
neling as described in Ref. [l^ due to the unequilibrated 
phonon effect. We ascribe this to the EPC induced selec- 
tive cascades of single-electron transitions between dif- 
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ferent charge states. In addition, we also analyze the 
current noise properties of this system. We find that 
the occurrence of weak and strong NDCs is accompanied 
by corresponding weak and strong enhancements of zero- 
frequency shot noise, respectively. 

To emphasize the external-voltage-driven unequili- 
brated phonon effect in electronic tunneling in a CNT, 
we consider the IVM in our model coupled to an exter- 
nal dissipative environment (a phonon bath),— and we 
incorporate the dissipation mechanism of the unequili- 
brated phonon into the ensuing rate equations on a mi- 
croscopic basis. For extremely strong dissipation, our 
results naturally reduce to those of equilibrated phonon- 
mediated tunneling. Therefore, our rate equations pro- 
vide a valuable theoretical framework for analysis of the 
mechanism by which finite phonon relaxation influences 
the electronic transport properties of CNT. 

The outline of the paper is as follows. In Sec. II, we de- 
scribe the model system that we study, and rewrite the 
model Hamiltonian in terms of an electron-phonon di- 
rect product (EPDP) state representation, which is suit- 
able for the ensuing theoretical derivation. In Sec. Ill, 
we derive a set of rate equations using a generic quan- 
tum Langevin equation approach with a Markovian ap- 
proximation, which facilitates investigation of uneqili- 
brated phonon and phonon dissipation effects on phonon- 
mediated electronic tunneling. In this section, we also 
derive the current formula using linear-response theory. 
In Sec. IV, we describe MacDonald's formula for calcu- 
lating zero-frequency shot noise by rewriting the ensuing 
rate equations in a number-resolved form. Then, we in- 
vestigate in detail the vibration-mediated transport and 
shot noise properties of a molecular QD with two MOs 
in Sec. V. Finally, a brief summary is given in Sec. VI. 



II. MODEL HAMILTONIAN 

In this paper, we consider a generic model for a molec- 
ular QD (in particular the suspended CNT) with two 
spinless levels, one as the highest-occupied MO (HOMO) 
El and the other as the lowest-unoccupied MO (LUMO) 
£2, coupled to two electrodes left (L) and right (R), and 
also linearly coupled to an IVM of the molecule having 
frequency with respective coupling strengths Ai and 
A2. We suppose that this single phonon mode is coupled 
to a dissipative environment represented as a set of inde- 
pendent harmonic oscillators (phonon-bath) . The model 
Hamiltonian is 

H = Hieads + Hmol + Hb + Hi, (la) 

with 

77, k 

Hmol = gjCjC, +Unin2 
i=i,2 



+uJoa'^a+ ^ Xjc'jCj{a'^ + a), (Ic) 

Hb = Y.^'v^lK^ (Id) 

p 

Hi = Ht + H,b, (le) 

Ht = ^(K„4c,+H.c.), (If) 

HvB = {a^ + a)^Kp{hl + bp), (Ig) 
p 

where c|^]^ (cjjk) is the creation (annihilation) operator of 
an electron with momentum k, and energy e^k in lead 
?7 (77 = L,R), and c] {cj) is the corresponding oper- 
ator for a spinless electron in the jth level of the QD 
(j = 1,2). U denotes interdot Coulomb interaction and 
Uj — CjCj is the electron number operator in level j. 
(a) and 6j (bp) are phonon creation (annihilation) op- 
erators for the IVM and phonon-bath (energy quanta 
loq, ujp), respectively. Xj represents the coupling con- 
stant between electron in dot j and the IVM; Kip IS the 
coupling strength between the IVM and phonon bath; 
V,fj describes the tunnel-coupling between electron level 
j and lead 77. Here, we denote the density of states of 
the phonon bath with respect to frequency cOp by D(ujp), 
and define the corresponding spectral density Jsi^p) as 

Jb(wp) = KlD{ujp). (2) 

In the literature, the following form of the spectral den- 
sity is usually considered j2S 

Jsiujp) = Koujp9{ujp), (3) 

in which 9{x) is the Heaviside step function. The bath is 
said to be Ohmic, sub-Ohmic, and super-Ohmic if a = 
1, a < 1, and a > 1, respectively. We use units with 
H — = e = 1 throughout the paper. In addition, it is 
worth noting that we do not consider the EPC-induced 
coupling between two MO's in the Hamiltonian, Eq. (fTc)l . 
because it is much weaker than the coupling of a single 
MO. 

It is well-known that, in the strong electron-phonon 
interaction problem, it is very convenient to introduce a 
standard canonical transformation^^ to the Hamiltonian 
Eq. (Hil), H = e^He~^, with S = {giiii + g2n2){a^ - a) 
(gj — Xj/ujo), which leads to a transformed Hamiltonian 

Hmol = "^ejcyj + Unin2 + uiQa^a, (4a) 
j 

Hi = Ht + Hvb, (4b) 

Ht = ^(T^,,4c,X,+H.c.), (4c) 

H^B = {a'' +a~2^gjnj)^Kp{bl + bp), (4d) 
j p 

with Sj = Ej - ^ and U = U - 2^^. Importantly, 
the X-operator describes the phonon renormalization of 
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dot-lead tunneling, 



„gj(a-a+) 



(4e) 



Obviously, the model described by the above Hamilto- 
nian Eq. (jlap involves a many-body problem with phonon 
generation and annihilation when an electron tunnels 
through the central region. Therefore, one can expand 
the electron states in the dot in terms of direct product 
states composed of single-electron states and n-phonon 
Fock states. In this two-level QD system, there are a 
total of four possible electronic states: (1) the two lev- 
els are both empty, |0)i|0)2, and its energy is zero; (2) 
the HOMO is singly occupied by an electron, |l)i|0)2, 
and its energy is ei; (3) the LUMO is singly occupied, 
|0)i|l)2, and its energy is £2; and (4) both orbitals are 
occupied, |l)i|l)2, and its energy is ei + 62 + U. Of 
course, if the interdot Coulomb repulsion is assumed to 
be infinite, the double-occupation is prohibited. For the 
sake of convenience, we assign these Dirac bracket struc- 
tures as dyadic kets (bras) )^ namely, the slave-boson kets 
4 = |0)i|0)2 «) 4 = |l)i|l)2 (E) |n), and pseudo- 
fermion kets fl = |l)i|0)2 <E) \n), /t„ = |0)i|l)2 ® \n). 
Correspondingly, the electron operator Cj and phonon op- 
erator a can be written in terms of such ket-bra dyadics 
as: 

00 

00 

a = Vn + l(et e„+-^ + ^ f^Jj^+i + dUn+iM 

with j 7^ j and Si(2) = ±1 (s2 = ^1 is due to anti- 
commutation relation of the Fermion operator). With 
these direct product states and dyadics considered as the 
basis, the density-matrix elements may be expressed as 
P^o = iPoo) = (4e„), p« = {p^/} = (4„/^.J, and = 

(/5Sd) = (4rf„>- 

The transformed Hamiltonian can be replaced by the 
following form in the auxiliary particle representation: 



rnol 



E 



^(ej +nwo)/j„/j, 



+ (ei + S2 + U + nLL>o)dld„], 



(7a) 



(7b) 



HyB = Kp{bl + 6„p) 



E Vn-I- l(ete„+i 



J2 fjnfjn+l + didn+l + H.C.) 



-2E5.(4/.n+4^^n) 



(7c) 



Based on this transformed Hamiltonian, we derive a rate 
equation for description of the dynamics of the reduced 



density matrix of the combined electron and IVM system 
in the sequential tunneling regime. 



III. QUANTUM LANGEVIN EQUATION 
APPROACH AND RATE EQUATIONS 

In this analysis, we employ a generic quantum 
Langevin equation approach )^^'^'^i^^i^^i^''''^^i^^ starting 
from the Heisenberg equations of motion (EOMs) for the 
density-matrix operators /5qq, {j = {1,2}), and p^^- 



ri^]s.,j,rn 



X (\/n + le,\e„+i + \/^e],e„_i - H.c), (8a) 

-v;,2rfL/i"^2^i™)sk] + E + ^-p) 

P 

x(Vn+ l//„/i„+i + ^/^flJm^i - H.c), 



(8b) 



«P22 



^Pdd 



+ Vjj2elnf2nX2,mn) + {Vill fti^mX^ ^^^ 
+ Vr,2dlf2nXl 

] + E'^p^^p + ^-p) 

P 

x(Vn+ l/]„/2n+i + Vnflj2n-i ~ H.c), 

(8c) 

[did^.H]^ = Y [cUV,,fLdnX,.,^r. 
r;,k,m 

-VTj2flmdnX2,mn) - (Kjl dJ,/2m-'^J,™„ 

-Vr,2diflmXl^Jc^]4 + Y '^P(K + ^-P) 

P 

X {Vn + Idid^^i + Vn4cf„-i - H.c), (8d) 



where 



X. 
X 



j.nm 

t 

j.nm 



(n|X,(t)|m), 
{m\X]{t)\n). 



(8e) 
(8f) 



These matrix elements can be calculated ass^i 



Xj,nm — Xj^^jj^ 



(9) 



where L™(a;) is the generalized Laguerre polynomial. 
The rate equations are obtained by taking statistical 
expectation values of the EOMs, Eqs. ([5a|) - (j8dp . which 
clearly involve the statistical averaging of products of 
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one reservoir (phonon-bath) variable and one device vari- 
able, such as X = (Lk.j, ^jji^k^n/ini^i,"™) and J = 

(Ep ^p{^\ + fe-p)\/^^(4e„+i - H.c.)). 

To determine the products, X and i/, we proceed by 
deriving EOMs for the system, phonon-bath, and reser- 
voir operators, Fi,„™ = ej^/^^, p^^^^ = 4 e«+i> and 



-[-Fl,nm, -f^r]- + [-Fl,nm, -H^ub]-, (10a) 

«Poo = [ej,e„+i,i7|_ = cjoPoo + -^tJ- 
+ ^ (10b) 

i5p = = t^p6p + [5p,i/„B]_, (10c) 

«c^k = Kk,^]- = e,7kc,,k + [c,7k,-^T]-- (lOd) 

The EOM for c^j^. is easily obtained by Hermitian con- 
jugation of the equations for c^\s_. Formally integrating 
these equations, pOap - (|10dp . from initial time to i we 
obtain 

Fi,n™(t) = e-^[^~i+(™-")"''l*^^i,„™(0) 

[ dt'e-'[^^+(™-")'^°l^[Fi,„„(t'),HT(t')]- 



Jo 



-ilei+{m-n)uJo] 



(11a) 



PoiT^'it) = e--*pj;6"+i(0) 



-z / dt'e-'-«^[pZ^\t'),HMt')]-, 



e-*'^''*6p(0) - i / dt'e-*"!- 

_ "'0 



bp{t) 



cr,k{t) = e-"'"'*c^k(0)-z j dt'e^^'^-^ 
xKk(0,^T(i')]-, 



(lib) 
(11c) 
(lid) 



with T ~ t~t' . In the absence of tunnel- coupling. Hi — > 
0, we have 

Knmit) = e-'[^^+(™-")-ol-F°_(t'), (12a) 
Poo^'W = e--°^pX+i(t'), (12b) 
fep(i) = e-^-^^bl{t'), (12c) 
c°k(0 = e-'"'-c°k(t')- (12d) 



A standard assumption in the derivation of a quan- 
tum Langevin equation is that the time scale of decay 
processes is much slower than that of free evolution, 
which is reasonable in the weak-tunneling approxima- 
tion. This bespeaks a dichotomy of time-developments of 
the involved operators into a rapidly- varying (free) part 



and a slowly- varying (dissipative) part. Focusing atten- 
tion on the slowly-varying decay processes, and noting 
that the infinitude of macroscopic bath variables barely 
senses reaction from weak interaction with the QD, it is 
appropriate to substitute the time-dependent decoupled 
reservoir, phonon-bath, and QD operators of Eqs. (|12ap - 
(|12dp into the integrals on the right of Eqs. I|llap - (|lldp . 
This yields approximate results for the reservoir opera- 
tors asj^s^^i^s 



c,k(i)=c°k(t)+c:^(i) 



(13a) 



with 



'^?7k( 



{t)^-z I dT[c;^{t),H^{t')]^, (13b) 



where H^^^gs^ is composed of the operators in Ht(vB) 
which are replaced by their decoupled counterparts (in- 
teraction picture). In fact, this is just the operator for- 
mulation of linear response theory. Similarly, the ap- 
proximate results for the QD and phonon-bath arc also 
divided into two parts: 



Knra{t)~^ / dr [F° (t) , i?^ (i')] - 

Jo 



dT[Fl^^it),H:,s{t')]-, 

/o 

Poo (*) = POO.O (t) + PQO,rT W + PoirvBW 



(14a) 



= Poo7\t)-^ / rfr[p^o7^(0,H^(t')]- 

"'O 

fdr[p-C\t).HUt')]~, (14b) 

"'O 

bp{t) = b;{t) + b;^y) 

= b;{t)-t f dT[bl{t),H:s{t')]-. (14c) 

Note that we use the super(sub)scripts rT and rvB de- 
note the reactions from tunnel-coupling and environmen- 
tal dissipation, respectively. 

Employing the approximate solutions of Eqs. (|13ap . 
(fT3b| and (fMaP . we can evaluate T as 

r;,k 



77, k 



(15) 



The statistical averages involved here can be taken sep- 
arately in regard to the electron ensembles of the reser- 
voirs (many degrees of freedom) and in regard to the 
few degrees of freedom of the QD EPDP states. Ac- 
cordingly, the statistical average of the product of one 
reservoir operator and one system operator factorizes in 
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the averaging procedure. Therefore, the statistical aver- 
age of c°|^(t)F°„„(t) vanishes. Moreover, in the sequen- 
tial picture of resonant tunneling, the tunneling rates 
and current are proportional to second-order tunnel- 
coupling matrix elements. We thus neglect the term 
'^r)k(^)-^fnm(0 it is proportioual to the third-order 
tunnel-coupling matrix element, 0(V^j). The other in- 
teraction terms arise from tunneling reaction, and are of 
second-order of V,^j . After some lengthy but straightfor- 
ward algebraic calculations, we obtain 



Poo 



"[1 fvi^)]pTl}-^l,nm-^l,nm^ 



in which 



V3 I 



(16a) 



(16b) 



denotes the tunneling strength between the molecular or- 
bital j and lead rj (g^ is the density of states of lead 77), 
/,;(e) is the Fermi-distribution function of lead 77 with 
temperature T, and 



Xj^nm = {n\Xj{t')\m), 
XL^ = {MX]{t')\n). 



Considering a{t) — e 
X 



'^a{t'), we have 



X] 



-z(m— n)a;oT 



X 



(16c) 
(16d) 



(16e) 
(16f) 



In the derivation of Eq. (|16ap . we assume that states 
with different phonon-numbers are completely decoher- 
ent, Poo(dd) = Poo{dd)^"m and p]p = p"/™™, owing to big 
energy difference between the two direct product states 
|j, n) and |j, m) if n ^ m. Moreover, a Markovian ap- 
proximation will be adopted by making the replacement 



dr 



dr, 



(17) 



in the statistical averaging, since we are interested in the 
long time scale behavior of these density matrix elements. 

Furthermore, J can be evaluated using Eqs. (|14b|) . 
pic)) with bp{t) = e-''^'-^bp{t') as 



J 



]KpV^([{b;Ht) + b;{t)) 



p 

ApICM 

Hb7''\t) + bl^^itWoC'it) - Poo,o"m) 
] + 1)Poo[("bH) + l)^(wp + Wo) 



P00,rvB\^)) 



+*E'*p("' + '^)Poo^Wb{uJp)5{uJp + ujq) 
p 

+ {nB{uOp) + l)5{uOp - uoq)], 



(18) 



in which nsi^Jp) = (e^J'/-^ — 1)~^ is the Bose-distribution 
function. With the assumption that the phonon-bath is 
composed of infinitely many harmonic oscillators having 
a wide and continuous spectral density, we can make the 
replacement, in the wide-band limit. 



dujpKlD{L,p){---), (19) 



yielding 



J = -in7pns(cjo)(n+l)/9oo+in7p(ns(wo) + l)("-+l)Poo"\ 

(20) 

with TUp — kqujq being constant due to Eq. ([3]). 

Following the same calculational scheme indicated 
above, we evaluated the other statistical expectation val- 
ues involved in the EOMs, Eqs. ((5a|) - ((8d)) . Finally, we ob- 
tained the following rate equations in terms of the direct 
product state representation of the density-matrix for the 
description of sequential tunneling through a molecular 
QD, accounting for the unequilibrated phonon effect and 
its modified tunneling rates, as well as IVM dissipation 
to the phonon environment: 



Poo = 



^Pll 



^P22 - {"^l 



^t,nm)Poo\ 



Pll 



P22 



E 



<^n + )P00 + ^„+lP00 + ^rt-lPoO '(21a) 
\j^l,mnPO0 ^ (^l,mn + ^2,nm)Pll + ^2,nmPlld\ 

771 

-(tn,+ + m-)p'i^ + m-^,p^+' + m+_,p^^\{21c) 

Pdd — E [^JnxriPll + ^t,mnP22 ~ (^l,7nn + ^2,mn)Pdd\ 
m 

-(tn+ -t- tn-)pS, -I- + w+_,p'i^\{21d) 

with the normalization relation J^niPoo + Pii + P22 + 
Pdd) = 1- The electronic tunneling rates are defined as 

Kj.nm = E ^VjlimfvC^J + {m - n)uJo), 

n 

(22a) 



E 

n 



E 



, F". 

/ J 7] J. Tim 



E ^viliml'^ - M^j + {m- n)cJo)], 
n 

yf+. 

/ J 7] J, Tim 



E ^vjllmfniU + + (»7i - n)uQ), 
n 

Vf- =Vf y 

/ y rij.nm / j VJ inm 

ri _ 
X[l - fr,{U +Ij + (TO-n)wo)], 



(22b) 



(22c) 



(22d) 



with the FC factor {p = min{m, n} and q = max{m, n}, 
denoting the smaller and larger of the quantities m and 



6 



n, respectively) 



him 



2 



„2 2|m-n|P! 



[Ll7-"l(g|)]^ (22e) 



+ (-'^^Ll,nm ^L2.mn)P22 LI ,mn ^ L2 .mn) Pddl ■ 

(26) 



describing the modification of tunnel- coupling due to 
phonon generation and emission during the electron tun- 
neling processes. This FC factor is symmetric, 7j'„„ — 
7^„, and satisfies the sum rules ^„ 7^^ = 7«m = 1- 
Obviously, these rates have specific physical meanings: 
r^nm (rj^ram) describes the tunneling rate of an elec- 
tron entering (leaving) level j with null occupancy of 
level J, together with the transition of vibrational quanta, 
n —f m [m n); while r^„„ (^~jnm) describes the tun- 
neling rate of an electron entering (leaving) level j with 
level j occupied, together with the corresponding transi- 
tion of the IVM state. 

The transition rates of phonon number states are 



„nB{uJo){n + 1), 
Wp(nB(wo) + l)n. 



(23a) 
(23b) 



which indicates that the state of the IVM changes from 
n to n -|- 1 (n — 1) by absorbing (emitting) a phonon 
from (to) the phonon-bath without change of the elec- 
tronic state. Note that + w~ defines the relaxation 
rate of the number state n of the IVM due to dissipa- 
tive coupling to the environment. Moreover, it should be 
pointed out that Wp — denotes no dissipation of the 
IVM to the environment, signifying the maximal unequi- 
librated phonon effect in tunneling; while Wp = 00 de- 
notes extremely strong dissipation of the IVM, so that it 
is always functions as an equilibrated state during each 
tunneling process, i.e., the excited phonon relaxes very 
quickly due to strong dissipation, before the next elec- 
tronic tunneling event takes place. Obviously, strong dis- 
sipation, Wp — 00, forces the probability distributions on 
the right-hand side of Eqs. (|21a[) - (j21dp to have the forms, 
Poo = Poo-P", P'jj = PjjP". and pJJ^ = RddP", with a 
thermal phonon distribution P" = e"""o/^(l - e^'^"/^) 
(assuming the phonon bath to have the same tempera- 
ture as the electrodes). 

The tunneling current operator through the molecular 
QD is defined as the time rate of change of the charge 
density, A^^ = I]k"],k"7,k' ™ I'^ad 77: 



(24) 



Employing linear-response theory in the interaction 
picture f2i we have 



I={Jl) 



dt'{[jL{t),m{f)]-)- (25) 



Following the same procedures indicated above ;^^'^^ we 
obtain the current formula through the left lead in terms 
of the density matrix elements of direct product states as 



L2,nm 



IV. MACDONALD'S FORMULA FOR 
ZERO-FREQUENCY SHOT NOISE 

In this section, we discuss the zero-frequency current 
noise of vibration-mediated sequential tunneling through 
a molecular QD involving an unequilibrated phonon. For 
this purpose, we employ MacDonald's formula for shot 
noise"^*^ based on a number-resolved version of the rate 
equations describing the number of completed tunneling 
eventsjSA This can be derived straightforwardly from the 
established QREs, Eqs. (PTa|) -(|21d | ). We introduce the 

two-terminal number-resolved density matrices p"j''' (t), 
representing the probability that the system is in the 
electronic state |j) (j = {0, 1,2, d}) with n vibrational 
quanta at time t together with 1(1') electrons occupying 
the left (right) lead due to tunneling events. Obviously, 

Pjj^^'^ ~ ^iv pfj '^^ ^'^^ resulting two-terminal 
number-resolved QREs for the case of an unequilibrated 
phonon are: 
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The current fiowing through the system can be evalu- 
ated by the time rates of change of electron numbers in 
lead rj as 



1,1' 



(28) 
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where 



(29) 

is the total probabiHty of transferring electrons 
into the left(right) lead by time t and Iri = l{l') with 
ry = L{R). It is readily verified that the current obtained 
from Eq. (f28|) by means of the number-resolved QREs, 
Eqs. (P7a|l - ((27d)) . is exactly the same as that obtained 
from Eq. (|26p. The zero- frequency shot noise with re- 
spect to lead 77 is similarly defined in terms of P^-'' •'(t) 
as well>2ai30^ 



^,(0)^2- 



Ll' 



(30) 



To evaluate S',j(0), we define an auxiliary function 



G];{t) as 



(31) 



Ll' 



whose equations of motion can be readily deduced em- 
ploying the number-resolved QREs, Eqs (j27ap - (j27d|l . in 
matrix form: G''(i) = 7W,,G''(t) + Gr^pit) with G''(t) = 
(GjJ(,,G^i,G^2:G2rf)'^ and p{t) = {poo, Pii, P22, PddV 



[here G^ = (G^, G^, • • • )^ and p,, = (pO^, p]^, • • • )^ 
Mr, and Qjj can be obtained easily from Eqs. (|27aP " (|27dp . 
Applying the Laplace transform to these equations yields 



G"(s) = (sI-M^r^GriPis), 



(32) 



where p{s) is readily obtained by applying the Laplace 
transform to its equations of motion with the initial con- 
dition p(0) = Pst [pst denotes the stationary solution of 
the QREs, Eqs (PTa|) - (|21d|) ]. Due to the inherent long- 
time stability of the physical system under consideration, 
all real parts of nonzero poles of p(s) and G^{s) are neg- 
ative definite. Consequently, the divergent terms aris- 
ing in the partial fraction expansions of p{s) and G^{s) 
as s — > entirely determine the large-t behavior of the 
auxiliary functions, i.e. the zero-frequency shot noise, 
Eq. jSOD. 

It is worth noting that (1) our two-terminal number- 
resolved QREs, Eqs. (|27aP " (|27dp . facilitate evaluation 
of the bias-voltage dependent zero-frequency shot noise 
for arbitrary interdot hopping; (2) our calculations yield 
Sl{0)=Sr{0). 



V. RESULTS AND DISCUSSION 

We now proceed with numerical calculations of the 
current / [Eq. (|26p ]. the zero- frequency current noise 
5(0) and the Fano factor F = S{0)/2I for the two-MO 
model in order to explain the particular experimental 



data in the current-voltage {I-V) characteristics recently 
reported for a suspended CNT.- 

For this purpose, we set the parameters in our calcula- 
tions as: cjo = 1 as the energy unit, Tjh/Tli = 10"^, 
rL2/rLi = 1, Tr2/Tli = 0.1, gi = 1, and ei = 0, 
£2 — 4.0wo- For simplicity, we fix the energy of the 
ground MO to be zero and ignore the nonzero bias- 
voltage- induced energy shift of the MO, which can be 
achieved using gate voltage in the experiment. Different 
from Ref. [TgI . we choose a large asymmetry in the tunnel- 
ing rates of the first MO, Tfn/TLi = 10"^, and an inter- 
mediate electron-phonon coupling strength, gi. In partic- 
ular, numerical fits of the experimentally measured data 
for the I-V curves show gi ~ 1 for long CNTsii These 
two parameters are believed to be necessary for the ap- 
pearance of NDC in combination with the unequilibrated 
phonon condition, for the regime in which the ground MO 
is predominant in tunneling (this will be shown below).— 
Moreover, we also choose a large asymmetry in the tun- 
neling rates of the second MO, Tii2/TL2 = 0.1, which is 
responsible for the appearance of strong NDC when the 
excited MO is predominant in tunneling as pointed out 
in Ref. 19. Furthermore, we examine the effect of EPC 
in the excited MO on the strong NDC. 

In our calculations, zup ~ denotes no dissipation of 
the IVM to the environment, corresponding to the max- 
imal uneqilibrated phonon effect in resonant tunneling; 
increasing dissipation strength, zup > 0, describes the ac- 
tion of the dissipative environment as it begins to relax 
the excited IVM towards an equilibrium phonon state 
(in the limit Wp = 00). This parameter allows one to 
examine the continuous variation of the effect of the dis- 
sipative environment on the excited IVM, which is help- 
ful in developing a deep understanding of the underlying 
properties of vibration-mediated resonant tunneling and 
its fluctuations in CNTs>i^ Throughout the paper we set 
the temperature as T = 0.02^0 and assume that the bias 
voltage V is applied symmetrically, fiL = —pR = V/2. 



A. Strong NDC in current-voltage characteristics 

Figure 1 exhibits the calculated results without envi- 
ronmental dissipation, vjp = 0, for the I-V character- 
istic and the occupation probability of the excited MO, 
P22, in the case of strong Coulomb blockade U = 00. 
If the excited MO is weakly coupled to the IVM, for 
instance g2 = 0.01 (solid lines) in Fig. 1, we find that 
(1) the ground MO plays a dominant role in tunneling if 
V < 222 (— 8); (2) a series of equally spaced steps in volt- 
age appears in the I-V curve demonstrating vibration- 
mediated transport behavior, associated with a gradually 
reduced step height due to the FC factor in this voltage 
range'r^iiSiiii'i^ (3) a weak NDC and current peak oc- 
cur at the onset of each step beginning from the second 
step, which are ascribed to the combination of asymmet- 
ric geometry, Tri/Tli = lO'^fi^ and the unequilibrated 
phonon effect (this will be shown below); and (4) a strong 
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NDC suddenly emerges when the excited MO becomes 
active in tunnehng &t V = 2?2, which helps to explain 
the experimental data [Fig. 3(a) in Ref. H]- As pointed 
out previously by M. Hettler et ai, who studied nonlinear 
transport through a molecular QD without IVM/^ this 
current decrease is a combined effect of strong Coulomb 
blockade and weak tunnel- coupling of the excited MO 
to the electrodes. It is for this reason that we selected 
r_R2/rLi = 0.1 in our numerical calculations. When the 
bias voltage increases to activate the excited MO, an elec- 
tron can occupy this MO but does not easily tunnel out 
due to the weaker escape rate Tii2, as shown by the solid 
line in Fig. 1(b), thus blocking occupation of the ground 
MO because of the strong Coulomb repulsion. As a re- 
sult, the tunneling current depends mainly on the con- 
tribution of the excited MO, leading to a suppression of 
current magnitude by about a factor r|j2- Relaxing ei- 
ther of these two conditions results in the elimination of 
the strong NDC. For example, the inset of Fig. 1 plots 
the corresponding results for the case without Coulomb 
interaction. Figure 2 shows the I-V curve with increasing 
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Figure 2: Current-voltage characteristics for several values of 
increasing tunneling rate r/jg/rLi with (a) g2 = 0.01 and (b) 
(?2 = 2.0 in the case of strong Coulomb blockade, U = oo, and 
unequilibrated phonon, ti7 = 0. 



of EPC, as shown in Fig. 1(b). In the following, we will 
offer a, Heener theoretical exnlanatinn of these results. 
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Figure 3: Electron-phonon joint occupation probabilities, poo 
(a,c) and P22 (b,d), vs. bias voltage relevant to Fig. 1 param- 
eters; for (a,b) (?2 = 2.0 and for (c,d) g2 = 5.0. 



Figure 1: (a) Calculated current /, and (b) occupation prob- 
ability P22 as functions of bias voltage in the case of strong 
Coulomb blockade, U = 00, and unequilibrated phonon, 
CC7 = for 172 = 0.01, 2.0, and 5.0. Inset; current / vs. voltage 
V in the case U — 0. 

Interestingly, if the excited MO is also appreciably 
coupled to the IVM, a stepdown behavior with equally 
spaced width in voltage is superposed on the overall de- 
crease of current, leading to a slowdown of the original 
rapid-reduction of current and more peaks of the NDC. 
More interestingly, an advancing and/or a postponing of 
the current decrease are observed, depending upon the 
relative strengths of the EPCs of the ground and ex- 
cited MOs. Intuitively, this advancing and postponing 
can be ascribed to the corresponding behaviors of the 
occupation probability of the second MO in the presence 



22 




Oo 



Figure 4: Allowed transitions for Fig. 3 (see text for details). 

For illustrative purposes, it is helpful to examine bias 
voltage dependences of the electron-phonon joint occu- 
pation probabilities (EPJOPs), Pqq and ^22- We plot the 
calculated results in Fig. 3 for two cases with 52 = 2.0 
(a,b) and 32 = 5.0 (c,d). It is unnecessary to plot the 
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EPJOPs for the ground MO, Pi^, because once an elec- 
tron enters into the first MO from the left electrode, it 
will escape very rapidly to the right electrode due to 
the strongly asymmetric configuration, Tjh/Tli = 10'^ 
(We only consider V > in this paper). Therefore, 
all are nearly zero. In contrast to this, once an 
electron is injected into the second MO, it is effectively 
trapped in this MO due to the suppressed tunnel-out rate 
^R2/^Li = 0.1 in the present model under investigation. 
With regard to this consideration, we first discuss the 
results for g2 = 2.0. Fig. 3(a) clearly shows that the 
EPDP state |0)i|0)2 'S) \n) (for notational convenience, 
we use 0„ to denote this EPDP state and 2„ to repre- 
sent |0)i|l)2 \n) below) is occupied, and contributes to 
current when the Fermi energy of the left lead is equal 
to the energy of the corresponding direct product state, 
PL — nujo, as illustrated in the schematic energy diagram. 
Fig. 4. Opening a new channel will cause a decrease of 
the occupation probabilities of previous direct product 
states. Intuitively, one might think that the EPDP state 
|0)i|l)2®|0)(= 2o) is unoccupied until the bias voltage in- 
creases to V = 2e2 = S.Owo- From Fig. 3(b), however we 
(surprisingly) observe that the EPDP state 2o is actually 
becoming occupied even at V = 4.0wo, which is half of 
the conventional resonant tunneling value. Moreover, the 
EPDP states 2i and 22 are both becoming occupied start- 
ing at V — G.OujQ, albeit their corresponding resonant 
values should traditionally be = lO.Oo^o and 12.0ijJo 
respectively. In addition, differing from the voltage de- 
pendence features of Pqq, the opening of new channels 
involving the excited MO does not cause a reduction of 
the occupation probabilities of previous channels. These 
peculiar properties can be understood qualitatively in 
terms of phonon-induced cascaded single-electron tran- 
sitions as illustrated in Fig. 4. As pointed out by M.R. 
Wegewijs, et al, arbitrarily high vibrational excitations 
can in principle be accessed via cascades of single-electron 
tunneling processes driven by a finite bias voltage, sir 
transitions between the EPDP states is related to t 
variation of electronic energy and the change of phone 
number states. For instance, if the Fermi energy of t 
left lead is located between the EPDP states O2 and 
with increasing bias voltage, i.e., A.Oujo < V < G.Oujq, t 
EPDP state O2 becomes occupied as shown in Fig. 3( 
Moreover, the single-electron transition, O2 2q, in 
cated by the arrow in Fig. 4 is also permitted becai 
the bias voltage V = i.Oivo provides sufficient enei 
to activate this transition, £2 + (0 — 2)wo = 2.0tJo- 
a result, albeit the Fermi energy of the left lead is r 
aligned with the energy £2 — 4.0a;o of the EPDP st; 
2o, this state also becomes occupied [Fig. 3(b)], whi 
precedes the conventional resonance value V = 8.0t 
Likewise, when the bias voltage increases to V > 6.0t 
the EPDP state O3 becomes occupied, and concomitant 
the transitions O3 — > 2o, O3 — > 2i, O3 22, and e-v 
Oi —^ 2q are also permitted with differing transition ra 
depending on the FC factors (denoted by the differe 
widths of the arrows in Fig. 4) via the cascade tran 



tion mechanism. Therefore, we observe from Fig. 3(b) 
that the states 2i and 22 both become occupied starting 
at V = 6.0liJq. In contrast to the situation in Ref. [la . 
the back-transition 2„ — > Om is prohibited in the present 
model due to the above-mentioned trapping effect of the 
excited MO (stemming from the suppressed escape rate 
rfl2)- Consequently, we find an accumulated increase of 
occupation probabilities of the EPDP states with low vi- 
brational excitations up to a threshold value of bias volt- 
age, in which a considerably large number of channels are 
stimulated and become active. In sum, the EPC-induced 
unidirectional cascaded transitions are responsible for the 
advanced appearance of NDC at lower bias voltages than 
one might otherwise expect. 

We now turn to examine the mechanism of the post- 
poned appearance of NDC in the case of (72 = 5.0 ^ gi. 
Figures 3(c) and (d) show that the occupation probabili- 
ties of Pqq have a similar bias voltage dependence to that 
of the case of 32 — 2.0 ~ gi, but the situation is con- 
siderably different for ^32- obviously, P22 is nearly zero 
until the bias voltage increases up to IO.Owq: which is 
even higher than the resonance value, V = Te^ = S.Oa-io. 
Therefore, it is interesting to explore why the EPC- 
induced cascade mechanism of single-electron transitions 
does not work in this situation. Albeit the cascade mech- 
anism for the transition, e.g. O2 — > 2o, satisfies the usual 
resonance condition from the energetic point of view, the 
actual occurrence of this transition still depends on the 
relevant transition rate determined by the correspond- 
ing FC factor, 733. In Fig. 5, we show the FC factors 
from n,TO = to 10 for g2 = 2.0 and 5.0, respec- 
tively. Clearly, a nearly vanishing FC factor occurs 
if 52 — 5.0. The strong EPC strength even blocks the 
conventional resonant transition O4 2q. Only when a 
large number of transitions are opened by a sufficiently 
high bias voltage (e.g. V = IO.Owq here), can a signifi- 
cant occupation of P22 be accumulated, thus leading to 




02468 10 02468 10 



Figure 5: Tlie FC factors, 7^^ {n,m = 0-10), for (a) 32 = 2.0 
and (b) (?2 ~ 5.0. 

We now further examine the disappearance of NDC 
shown in Fig. 2(b). With gradual relaxation of the trap- 
ping effect of the excited MO (increasing Tii2), the unidi- 
rectional cascaded transitions becomes bidirectional, i.e., 
not only can the transitions 0„ — s- 2^ occur, but also 
2n Om occurs as well, provided that the energy conser- 
vation condition is satisfied and the FC factor permits. 
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Therefore, the inverse transition, 2q O2, reduces the 
occupation of the state 2q (not shown here) , and finaUy 
reduces the NDC. 



B. Effects of dissipation to environment and 
super-Poissonian current noise 

We now discuss the effects of environmental dissipation 
on the current and zero-frequency shot noise. Figures 6 
and 7 exhibit the dissipation dependences of the current / 
and Fano factor F = S{Q)/2I as functions of bias voltage 
for systems having 32 — 0.01 and 2.0, respectively. 
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Figure 6: Calculated current, / (a) and Fano factor, F — 
S{0)/2I (b) vs. bias voltage for the case of g2 — 0.01 and 
various environmental dissipation rates rop. Other parameters 
are the same as in Fig. 1. 
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Figure 7: The same figure as Fig. 6 except that 172 = 2.0. 

In the region where the ground MO is dominant in 
tunneling, it is clearly observed that the weak NDC 
becomes ever weaker with a gradually increasing dissi- 
pation rate tu. The FV characteristic finally exhibits 
only positive differential conductance at w = 00, i.e.. 



the equilibrated phonon condition. Therefore, one can 
conclude that the observation of weak NDC in the re- 
cent transport measurements clearly indicates that the 
external voltage-driven unequilibrated phonon effect in 
a suspended CNT plays an essential role in determin- 
ing its underlying transport properties;^ Moreover, it is 
also clear that, due to the unequilibrated phonon effect, 
the current noise shows a weak super-Poissonian char- 
acteristic {F > 1) associated closely with the appear- 
ance of NDC. Furthermore, the environmental dissipa- 
tion suppresses the Fano factor and for certain values 
of the dissipation rate, tu, the shot noise shows weak 
sub-Poissonian behavior {F < 1), but it becomes Pois- 
sonian {F = 1) at the equilibrated phonon condition. 
This is just the traditional value of the Fano factor, 
F = (rii + r|i)/(rii -t- Tri^, for the extremely asym- 
metrical configuration F/ji/Fli = 10'^. 

From Figs. 6 and 7, we find strong super-Poissonian 
shot noise as a companion to the strong NDC when the 
second MO starts to contribute to the current, irrespec- 
tive of the environmental dissipation. When the second 
MO is also coupled to the IVM, the environmental dissi- 
pation strongly influences the current and noise as shown 
in Fig. 7. In particular, the threshold value of bias volt- 
age for strong NDC and strong enhancement of shot noise 
is exactly the traditional resonant tunneling value of the 
second MO. 



VI. CONCLUSIONS 

In summary, we have fully analyzed the external-bias- 
voltage-driven nonequilibrated phonon effect on nonlin- 
ear tunneling through a suspended CNT in the sequential 
tunneling regime. In order to qualitatively address the re- 
cent experimental results, we have modeled the CNT as a 
molecular QD having two electronic MOs with asymmet- 
ric tunnel-coupling rates to the left and right electrodes 
and strong interaction with an IVM. To study the role of 
dissipation of unequilibrated phonons, we have assumed 
further that the molecular IVM is also weakly coupled to 
a phonon bath "environment" . To carry out this anal- 
ysis, we established generic rate equations in terms of 
the EPDP state and auxiliary-particle representation for 
the description of vibration-mediated resonant tunneling 
employing a microscopic quantum Langevin equation ap- 
proach in the limit of weak tunneling and weak dissipa- 
tion. 

Employing the ensuing rate equations derived here, 
we systematically analyzed vibration-mediated resonant 
tunneling in the present model, obtaining the FV char- 
acteristics, zero-frequency current noise, and the effects 
of environmental dissipation on the role of the unequi- 
librated phonons. Our numerical analysis shows that 
in the voltage region where the ground orbital is dom- 
inant in tunneling, the combined effect of unequilibrated 
phonons and asymmetric tunnel- couplings is responsible 
for weak peaklike structures in the FV curve at the on- 
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set of each phonon step with a weak NDC and corre- 
spondingly weak super-Poissonian noise. Furthermore 
we found that this peakhke structure could be gradually 
diminished by environmental dissipation of the uneqili- 
brated IVM and become completely devoid of NDC at the 
equilibrated phonon condition. Accordingly, the usual 
current noise for an asymmetric single-electron tunneling 
device is predicted at the equilibrated phonon condition, 
-F = 1; however, for a finite dissipation rate suppressed 
noise, F < 1, may be observed under certain conditions. 

More interestingly, we have also discussed the trans- 
port properties in detail in the second-orbital dominated 
region. We found that the interplay of strong Coulomb 
interaction between the two MOs and strong asymmetry 
of tunnel-coupling leads to very strong NDC and a cor- 
respondingly strongly enhanced shot noise, irrespective 
of the EPC. However, the bias voltage value for the on- 
set of NDC and super-Poissonian shot noise is intimately 
dependent on the EPC strength of the second MO, i.e., 
this voltage value can be smaller or larger than the tradi- 



tional resonant tunneling value for the second MO. Our 
discussion concluded that this feature stems from the 
EPC-induced selective cascades of single-electron tran- 
sitions with FC-factor-modified rates in a unidirectional 
way due to the asymmetric tunnel- coupling. Environ- 
mental dissipation forces this value to tend to the tradi- 
tional resonance point. 
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